CNV-association meta-analysis in 191,161 European adults reveals new loci associated with anthropometric traits

Funding Information: This research has been conducted using the UK Biobank Resource. This research has been conducted using the Danish National Biobank resource. The authors are grateful to the Raine Study participants and their families, and to the Raine Study research staff for cohort co-ordination and data collection. QIMR is grateful to the twins and their families for their generous participation in these studies. We would like to thank staff at the Queensland Institute of Medical Research: Anjali Henders, Dixie Statham, Lisa Bowdler, Ann Eldridge, and Marlene Grace for sample collection, processing and genotyping, Scott Gordon, Brian McEvoy, Belinda Cornes and Beben Benyamin for data QC and preparation, and David Smyth and Harry Beeby for IT support. HBCS Acknowledgements: We thank all study participants as well as everybody involved in the Helsinki Birth Cohort Study. Helsinki Birth Cohort Study has been supported by grants from the Academy of Finland, the Finnish Diabetes Research Society, Folkhälsan Research Foundation, Novo Nordisk Foundation, Finska Läkaresällskapet, Juho Vainio Foundation, Signe and Ane Gyllenberg Foundation, University of Helsinki, Ministry of Education, Ahokas Foundation, Emil Aaltonen Foundation. Finrisk study is grateful for the THL DNA laboratory for its skillful work to produce the DNA samples used in this study and thanks the Sanger Institute and FIMM genotyping facilities for genotyping the samples. We thank the MOLGENIS team and Genomics Coordination Center of the University Medical Center Groningen for software development and data management, in particular Marieke Bijlsma and Edith Adriaanse. This work was supported by the Leenards Foundation (to Z.K.), the Swiss National Science Foundation (31003A_169929 to Z.K., Sinergia grant CRSII33-133044 to AR), Simons Foundation (SFARI274424 to AR) and SystemsX.ch (51RTP0_151019 to Z.K.). A.R.W., H.Y. and T.M.F. are supported by the European Research Council grant: 323195:SZ-245. M.A.T., M.N.W. and An.M. are supported by the Wellcome Trust Institutional Strategic Support Award (WT097835MF). For full funding information of all participating cohorts see Supplementary Note 2. Publisher Copyright: © 2017 The Author(s).There are few examples of robust associations between rare copy number variants (CNVs) and complex continuous human traits. Here we present a large-scale CNV association meta-analysis on anthropometric traits in up to 191,161 adult samples from 26 cohorts. The study reveals five CNV associations at 1q21.1, 3q29, 7q11.23, 11p14.2, and 18q21.32 and confirms two known loci at 16p11.2 and 22q11.21, implicating at least one anthropometric trait. The discovered CNVs are recurrent and rare (0.01-0.2%), with large effects on height (> 2.4 cm), weight ( 5 kg), and body mass index (BMI) (> 3.5 kg/m(2)). Burden analysis shows a 0.41 cm decrease in height, a 0.003 increase in waist-to-hip ratio and increase in BMI by 0.14 kg/m2 for each Mb of total deletion burden (P = 2.5 x 10(-10), 6.0 x 10(-5), and 2.9 x 10(-3)). Our study provides evidence that the same genes (e.g., MC4R, FIBIN, and FMO5) harbor both common and rare variants affecting body size and that anthropometric traits share genetic loci with developmental and psychiatric disorders.Peer reviewe

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2,907 cases with AN from 14 countries (15 sites) and 14,860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery datasets. Seventy-six (72 independent) SNPs were taken forward for in silico (two datasets) or de novo (13 datasets) replication genotyping in 2,677 independent AN cases and 8,629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication datasets comprised 5,551 AN cases and 21,080 controls. AN subtype analyses (1,606 AN restricting; 1,445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01×10−7) in SOX2OT and rs17030795 (P=5.84×10−6) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76×10−6) between CUL3 and FAM124B and rs1886797 (P=8.05×10−6) near SPATA13. Comparing discovery to replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P= 4×10−6), strongly suggesting that true findings exist but that our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field

A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling

J. Kaprio, A. Palotie, A. Raevuori-Helkamaa ja S. Ripatti ovat työryhmän Eating Disorders Working Group of the Psychiatric Genomics Consortium jäseniä. Erratum in: Sci Rep. 2017 Aug 21;7(1):8379, doi: 10.1038/s41598-017-06409-3We conducted a genome-wide association study (GWAS) of anorexia nervosa (AN) using a stringently defined phenotype. Analysis of phenotypic variability led to the identification of a specific genetic risk factor that approached genome-wide significance (rs929626 in EBF1 (Early B-Cell Factor 1); P = 2.04 x 10(-7); OR = 0.7; 95% confidence interval (CI) = 0.61-0.8) with independent replication (P = 0.04), suggesting a variant-mediated dysregulation of leptin signaling may play a role in AN. Multiple SNPs in LD with the variant support the nominal association. This demonstrates that although the clinical and etiologic heterogeneity of AN is universally recognized, further careful sub-typing of cases may provide more precise genomic signals. In this study, through a refinement of the phenotype spectrum of AN, we present a replicable GWAS signal that is nominally associated with AN, highlighting a potentially important candidate locus for further investigation.Peer reviewe

Defining the Effect of the 16p11.2 Duplication on Cognition, Behavior, and Medical Comorbidities

IMPORTANCE: The 16p11.2 BP4-BP5 duplication is the copy number variant most frequently associated with autism spectrum disorder (ASD), schizophrenia, and comorbidities such as decreased body mass index (BMI). OBJECTIVES: To characterize the effects of the 16p11.2 duplication on cognitive, behavioral, medical, and anthropometric traits and to understand the specificity of these effects by systematically comparing results in duplication carriers and reciprocal deletion carriers, who are also at risk for ASD. DESIGN, SETTING, AND PARTICIPANTS: This international cohort study of 1006 study participants compared 270 duplication carriers with their 102 intrafamilial control individuals, 390 reciprocal deletion carriers, and 244 deletion controls from European and North American cohorts. Data were collected from August 1, 2010, to May 31, 2015 and analyzed from January 1 to August 14, 2015. Linear mixed models were used to estimate the effect of the duplication and deletion on clinical traits by comparison with noncarrier relatives. MAIN OUTCOMES AND MEASURES: Findings on the Full-Scale IQ (FSIQ), Nonverbal IQ, and Verbal IQ; the presence of ASD or other DSM-IV diagnoses; BMI; head circumference; and medical data. RESULTS: Among the 1006 study participants, the duplication was associated with a mean FSIQ score that was lower by 26.3 points between proband carriers and noncarrier relatives and a lower mean FSIQ score (16.2-11.4 points) in nonproband carriers. The mean overall effect of the deletion was similar (-22.1 points; P 100) compared with the deletion group (P < .001). Parental FSIQ predicted part of this variation (approximately 36.0% in hereditary probands). Although the frequency of ASD was similar in deletion and duplication proband carriers (16.0% and 20.0%, respectively), the FSIQ was significantly lower (by 26.3 points) in the duplication probands with ASD. There also were lower head circumference and BMI measurements among duplication carriers, which is consistent with the findings of previous studies. CONCLUSIONS AND RELEVANCE: The mean effect of the duplication on cognition is similar to that of the reciprocal deletion, but the variance in the duplication is significantly higher, with severe and mild subgroups not observed with the deletion. These results suggest that additional genetic and familial factors contribute to this variability. Additional studies will be necessary to characterize the predictors of cognitive deficits

Rare Genomic Structural Variants in Complex Disease: Lessons from the Replication of Associations with Obesity

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease

Chromosomal contacts connect loci associated with autism, BMI and head circumference phenotypes

Copy number variants (CNVs) are major contributors to genomic imbalance disorders. Phenotyping of 137 unrelated deletion and reciprocal duplication carriers of the distal 16p11.2 220 kb BP2-BP3 interval showed that these rearrangements are associated with autism spectrum disorders and mirror phenotypes of obesity/underweight and macrocephaly/microcephaly. Such phenotypes were previously associated with rearrangements of the non-overlapping proximal 16p11.2 600 kb BP4-BP5 interval. These two CNV-prone regions at 16p11.2 are reciprocally engaged in complex chromatin looping, as successfully confirmed by 4C-seq, fluorescence in situ hybridization and Hi-C, as well as coordinated expression and regulation of encompassed genes. We observed that genes differentially expressed in 16p11.2 BP4-BP5 CNV carriers are concomitantly modified in their chromatin interactions, suggesting that disruption of chromatin interplays could participate in the observed phenotypes. We also identified cis- and trans-acting chromatin contacts to other genomic regions previously associated with analogous phenotypes. For example, we uncovered that individuals with reciprocal rearrangements of the trans-contacted 2p15 locus similarly display mirror phenotypes on head circumference and weight. Our results indicate that chromosomal contacts’ maps could uncover functionally and clinically related genes

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 x 10(-7)) in SOX2OT and rs17030795 (P=5.84 x 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 x 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 x 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 x 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field